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ABSTRACT 

We present an analysis of four off-axis ROSAT PSPC observations of the Perseus 
cluster of galaxies (Abell 426). We detect the surface brightness profile to a radius of 
80 arcmin (~ Mpc) from the X-ray peak. The profile is measured in various 

sectors and in three different energy bands. Firstly, a colour analysis highlights a 
slight variation of iVn over the region, and cool components in the core and in the 
eastern sector. We apply the /3-model to the profiles from different sectors and present 
a solution to the, so-called, /3-problem. The residuals from an azimuthally-averaged 
profile highlight extended emission both in the East and in the West, with estimated 
luminosities of about 8 and f xlO'*'^ erg s~^, respectively. We fit several models to 
the surface brightness profile, including the one obtained from the Navarro, Frenk and 
White (1995) potential. We obtain the best fit with the gas distribution described by 
a power law in the inner, cooling region and a /?-model for the extended emission. 

Through the best-fit results and the constraints from the deprojection of the 
surface brightness profiles, we define the radius where the overdensity inside the cluster 
is 200 times the critical value, r2oo, at 2.7h^Q Mpc. Within 2.3h^Q Mpc (0.85r2oo), 
the total mass in the Perseus cluster is 1.2 x IO^^Mq and its gas fraction is about 30 
per cent. 

Key words: galaxies: individual: A426 - cooling flows - dark matter - X-ray: galaxies. 



1 INTRODUCTION 

The Perseus cluster (Abell 426) is a Bautz-Morgan type II- 
III cluster, of richness class 2 , with a cD galaxy (NGC 1275) 
at its centre (02000 3''19'"46.9", <52ooo 41°30'48"). The cluster 
is located at redshift 0.0178 (Fadda et al. 1996; 1 arcmin ~ 
30.1 kpc, using hereafter Ho = 50 km s^^ Mpc^^, qo — 0.5 
and A = 0). 

The Perseus cluster is the brightest cluster in the X- 
ray waveband, and so has been extensively studied since the 
UHURU mission (Gursky et al. 1971, Forman et al. 1972). 
The emission from its intracluster medium (ICM) is well 
described by a thermal gas model with a metallicity of 40 
per cent of the Solar value and a temperature of 6.3 keV (cf. 
Allen et al. 1992, using the GINGA satellite). The Spartan 1 
detector (Snyder et al. 1990; Kowalski et al. 1993) indicates 
an isothermal temperature profile, beyond the cooling flow 
region, out to 50 arcmin radius. 

However, more recent work shows that the ICM has a 
complex morphology in X-rays. The 505 s exposure from 
the ROSAT AW Sky Survey (Schwarz et al. 1992) has shown 
cooler emission (by a factor of 2.5-3 times lower with respect 
to the rest of the cluster) to the East, corresponding to a sus- 
pected merger (recognizable between 10 and 78 arcmin from 



NGC 1275). An analysis of the moments of the X-ray sur- 
face brightness distribution from the Einstein Observatory 
IPC observations provides additional evidence that the clus- 
ter is not yet in a relaxed state (Mohr, Fabricant & Geller 
1993). Furthermore, analysis of the ASCA CIS data (Ar- 
naud et al. 1994) indicates hotter regions to the North- West 
and South- West directions, and also a slight increases in the 
metal abundance. A negative linear gradient in the metallic- 
ity has also been observed in the Spartan 1 data (Kowalsky 
et al. 1993). Furthermore, a cooling flow with a deposition 
rate larger than 200 M0yr~^ is located around NGC 1275, 
which has an active, radio-loud nucleus (Fabian et al. 1981, 
Branduardi-Raymont et al. 1981, Primini et al. 1981, Roth- 
schild et al. 1981). 

In this paper, we report on the results of surveying the 
cluster out to radii of 80 arcmin (~ 2.4 Mpc) using a mo- 
saic of four off-axis, ROSAT Position Sensitive Proportional 
Counter (PSPC) observations. 

We present in section 2 the data and analysis of the 
surface brightness profiles using both a deprojection tech- 
nique and a detailed colour analysis. We also attempt to 
model the gas profile with a /3-model, and discuss the dif- 
ferent behaviour according to the regions and energy bands 
considered. 
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Figure 1. The unfolded spectrum is shown with its various 
additive components (up panel) and best-fit residuals (bottom 
panel). The spectrum has been made collecting the 0.2-2 keV 
photons from the central 20'— radius region in the West region. 
The adopted model is the sum of a Galactic emission modelled 
by a Galactic Raymond-Smith thermal plasma with Raymond- 
Smith cluster emission, with metallicity of 40 per cent of the 
Solar value, folded through a Galactic absorption of 12.6 xlO^'' 
atoms cm~^. The best-fit (x^ of 211 for 173 degrees of freedom) 
has parameters: fcTcai = 0.221°-°^ keV, /cT^iustcr = lO.iglj'j? 
keV. 
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Figure 3. An azimuthal plot of the surface brightness within 
sectors of 22?5, extracted around the X-ray peak in the on-axis 
image in steps of 4 arcmin for the region enclosed within < 
r < 24 arcmin. For comparison the values in ordinate for the 
azimuthal profiles extracted at r > 4' are multiplied by a factor 
of 3. Due north is defined as 0° and East at 90° . 
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Figure 2. A mosaic of the raw data of the 4 offset observations 
of the Perseus Cluster, selected in the hard energy band (52-201 
PI channels) and exposure-corrected, is shown here with an adap- 
tively smoothed contours. The eight contours are equally spaced 
in a logarithmic scale, starting from 1.5 xlO~^ cts arcmin"^, 
which is 3 times the averaged background of the four images, and 
increasing by a factor of 2 to 0.192 cts s~^ arcmin"^. 



Our purpose is to place strict constraints on the non- 
luminous matter in the cluster, under the assumption of 
hydrostatic equilibrium, as discussed in section 3. The main 
conclusions are summarized in section 4. 



2 X-RAY DATA AND ANALYSIS 

In Table 1, we summarise the ROSAT public archive obser- 
vations that we consider in our analysis. 

The Perseus cluster centre is located at Galactic coordi- 
nates {l,b: 150? 58, -13?26). Because of its low Galactic lati- 
tude. Galactic emission provides a dominant contribution at 
very low energy (~ 0.1 keV), making (i) the corresponding 
energy band of little use, and (ii) any attempt to constrain 
fits of the spectra extracted from the off-axis regions diffi- 
cult (for example, consider Fig. 0). 

Following Snowden et al. (1994), we correct the images 
by their exposure maps in the 0.5-2.0 keV band using the 
Interactive Data Language (IDL) ROSAT user-supplied li- 
braries. A mosaic of the four off-axis observations is shown 
in Fig. ^, together with adaptively smoothed contours of the 
exposure-corrected image. 

This figure clearly shows asymmetric emission along the 
East-West axis, which is also apparent in an azimuthal plot 
of the counts extracted at different radii from the central 
peak (cf. Fig. ||). The eastern excess (around position angle, 
p.a.= 90°) becomes significant at radii larger than 4 arcmin, 
whereas a larger contribution in emission along the West- 
Southwest direction (p. a. ~ 260°) is evident between and 
4 arcmin, coinciding with the position given by Schwarz et 
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Figure 5. The surface brightness profiles in bins of 4 arcmin ware extracted from the corresponding ofi-axis image and exposure-corrected 
for the counts in the 0.5-2.0 keV band. The dotted line indicates the background value as estimated in the 45°-sector opposite to the 
peak of the X-ray emission. The over-plotted points are the surface brightness profiles extracted, in bins of width of 3 arcmin, from the 

respective sectors of the on-axis image. In this case, the background (dashed line) has been estimated 40 arcmin from the X-ray peak in 
emission. The statistical error in each bin of the surface brightness profile is smaller than the dimension of the point symbol adopted. 



Table 1. ROSAT observation summary 



Seq. Id. Date Exp (s) 0:2000 <52000 Region N: 



ON-axis 

wp800186 1992 02 Feb 4787 3'U9™ 48?0 +41°30'36" Centre 14.9 

OFF-axis 

rp800035 1991 18 Aug 23669 3''18™ 28?7 -|-42°18'36" North 15.2 

rp800032 1991 02 Sep 29640 3''15'" 33!6 -|-41°16'48" West 12.6 

rp800034 1991 19 Aug 23348 3^21'" 09?6 -|-40°44'24" South 14.3 

rp800033 1991 09 Feb 14566 3''24™ 02?4 -|-41°46'48" East 17.0 

i This value (in units of 10'^'^ atoms cm~^) is the average of the HI column density, 
calculated by interpolating the {a, S) coordinates from the l°grid maps 
by Dickey & Lockman 1990. 
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Figure 6. Here all the background-subtracted ofT-axis surface brightness profiles, selected in the 0.5-2.0 keV band, are plotted together 
(using the corresponding symbol from Fig. S) with their average (Centre-Ojff profile; thick solid line). 



al. (1992; cf. their Fig. 2). The latter enhancement is also 
seen in a wavelet analysis of the on-axis ROSATPSPC image 
(cf. fig. 10 in Slezak, Durret & Gerbal 1994), and continues 
in the second bin (4-8 arcmin). This supports the idea of 
an asymmetric distribution in the gas, with a preferential 
orientation along the 'chain' of bright galaxies to the West, 
out to a radius of 10 arcmin (cf. the optical image of the 
central part of the Perseus cluster in Fig. and with the 
main excess above 10 arcmin to the East. Here, we note 
that further support for the idea of a merger proceeding 
along East-West axis is provided by the higher intracluster 
gas temperature observed by ASCA (Arnaud et al. 1994) in 
the North- West and South- West regions, where the effect of 
the shocks are expected to increase the cluster temperature 
(cf. the results of hydrodynamic simulations in Schindler & 
Miiller 1993 and Roettinger, Loken and Burns 1997). 

Moreover, an elliptical-fitting analysis with the TRAP 
routine ellipse shows the same behaviour of the p. a. shown 
in Fig. ^. We measure an ellipticity that increases from 
0.06 ± 0.04 at 4 arcmin to 0.22 ± 0.07 at 24 arcmin (be- 
yond this radius the error becomes larger than 50 per cent 
the value). The centroid oscillates around the X-ray peak 
by about 2 arcmin, first slightly to the West, then to the 
East, and finally Westward again (cf. Fig. Apart from 



the Perseus cluster (cf. also Mohr et al. 1993), such oscilla- 
tions are also seen in other clusters (eg. A2256, Mohr et al. 
1993; Centaurus cluster, Allen & Fabian 1994; A478, White 
et al. 1994; consider also the studies on the behavior of the 
X-ray isophotes for a sample of clusters studied by Buote & 
Canizares 1996). These are indicators of an on-going merger 
that causes a disturbance in the intracluster gasrl By study- 
ing the residuals of the image, once the best-fit elliptical 
isophotal model is subtracted, we find significant structure 
only if the centroids of the isophotes are fixed to the position 
of the X-ray peak. Allowing each isophotes centroid to vary, 
eliminates significant residuals. In the following analysis, we 
will compare our data against the fixed centroid model of 
the cluster. 

To extract the surface brightness profiles studied in the 
following analysis and presented in Figs. ^ and ^, each obser- 
vation has had the region of the PSPC support ribs masked 
out, as well as all detected point sources. Avoiding the cen- 
tral 1 arcmin where the non-thermal emission contributed 



* Also an optical analysis shows a strong dependence of the po- 
sition of the cluster center upon limiting magnitude and galaxy 
morphological types (Andreon 1994 and references therein) 
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Figure 4. The adaptively smoothed contours, with the same 
levels as in Fig. ^, are here plotted on the Digitized Sky Survey 
image of the central (30 arcmin)-^ area of the Perseus cluster. 
The chain of bright galaxies is evident westwards, in a direction 
of NGC 1272, located at about 5 arcmin from the cD, NGC 1275. 
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Figure 7. The eight contours of Fig. ^ are plotted on the 3° X 
'i° IRAS 100 micron map. The IRAS data are representative of 
the Galactic absorption and show, on the left, the plane of the 
Galaxy at latitude 6 ~ —12°. The (minimum, maximum) values 
in the IRAS map are ( 6.1, 19.2) xlO^ Jy sr"!. 



Figure 8. The Binney & Tremaine (BT) and Navarro-Frenk- 
White (NFW) dark matter profiles applied in the deprojection 
analysis are here compared for different input parameters, [a (km 
s~^), Tc (Mpc)]. Both of them are also compared with the King 
approximation to the true isothermal sphere (King 1962). All of 
them are normalized to the central value of the true isothermal 
sphere profile [po = 9cr2/(47rGr2) = 9.05 X 10"^^ g cm'^]. In- 
side the core radius, the NFW profile does not flatter like the 
BT profile. In the outer part of the region of interest (2-3 Mpc), 
agreement between the two profiles is obtained by increasing the 
velocity dispersion and the core radius (or scale radius) in the 
NFW profile. Here it is also shown the dependence on r~^'*, as 
recent N-body simulations of clusters of galaxies have obtained 
(cf. Thomas et al. 1998; the slope of this power law is the mean 
value estimated in the range 0.1-1 rigo; where rigo is the radius 
within that the average cluster density is 180 the critical one). 



by NGC 1275 could be significant, the surface brightness 
profile iias been extracted witii bins of width of 3 and 4 ar- 
cmin, for the on-axis and the ofF-axis images, respectively. 
These bins widths allow smoothing over the masked circular 
rib at about 20 arcmin for the on-axis images, and the Point- 
Spread-Function effect in the off-axis region of the ROSAT 
PSPC (cf. Fig. 5e in Hasinger et al. 1993). The profiles have 
been extracted from four 90° sectors, North, West, South 
and East, respectively, of the single on-axis image. Each off- 
axis observation provides an azimuthally averaged surface 
brightness profile. We have also considered an average of 
these profiles (called Centre-ojf). 

The background has been estimated in a 45° sector that 
avoids the ribs, in a region between 40 and 45 arcmin from 
the centre of the image: South-East in the on-axis image, 
and from the opposite side of the field, from the peak of 
the cluster, in the offset images. The X-ray peak has been 
estimated by applying on the correspondent image a median 
filter by a square moving box with width of 5 pixels (cf. 
Table 2). 

In Fig. P we show a comparison between the on-axis 
and the off-axis profiles. Due to the limited field-of-view of 
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Table 2. X-ray centres and their separation from NGC 1275 



Region a2000 52000 



NGC 1275 




48! 16 


+41 ° 30' 42'.' 1 




On-axis 




46! 89 


+41°30'48'.'5 


0;"33 


North 


3h;^gm 


46! 97 


+41°32'11'.'8 


li^52 


West 


19™ 


42! 65 


+41°30'23'.'5 


1™41 


South 


3h]^gm 


43! 51 


+41°29'59'.'l 


1™37 


East 


3h]^gm 


49! 31 


+41°31'37'.'9 


K"^97 



the on-£ixis image and the consequential diflerence in esti- 
mates of the background, we only consider the off-axis sur- 
face brightness profile in the following analysis. 

A comparison of the profiles shown in Fig. ^ indicates 
two enhancements with respect to the averaged value: (a) in 
the eastern sector out to 25 arcmin ; (b) beyond 30 arcmin 
in the West. Galactic absorption over the region (cf. Fig. ^ 
implies a discrepancy, in the ROSAT band flux, of up to 
17 per cent between East and West, for the most and least 
absorbed emission respectively. This compares with the dif- 
ference in the surface brightness values of a factor of about 
2 at l°radius from the X-ray centre. This underlines that 
the elliptical shape of the X-ray emission is not primarily 
due to the gradient in foreground absorption present over 
the region. 

A further discussion about the intrinsic reality of these 
variations is presented in sections 2.2 and 2.3. 

2.1 The deprojection analysis 

A deprojection analysis has been applied to the surface 
brightness profiles extracted from exposure-corrected im- 
age, with photons selected in the 0.5-2.0 keV band. (The 
original deprojection work by Fabian et al. 1981 was on the 
Perseus cluster; for recent details on the deprojection tech- 
nique, see White et al. 1997). Given the assumption that the 
observed projected cluster emission is due to the X-ray emit- 
ting gas in spherically symmetric shells, the count emissivity 
in each radial volume shell can be calculated analytically and 
compared with the predicted counts from a bremsstrahlung 
emission of an optically thin gas [described by a MEKAL 
model, based on the model calculations of Mewe and Kaas- 
tra (Kaastra 1992) with Fe L calculations by Liedahl (1995)], 
absorbed by intervening matter (A^h from Table 1) and con- 
volved with the response of the detector. For such an emis- 
sion model, the calculated fiux is proportional to the electron 
density, n^, and the intracluster medium temperature, Tgas, 
according to the relation: -rxifiT^^s) (where /(Tgas) ^ T^i^, 
for Tgas 3 X 10^ K) . Applying the perfect gas law and the 
equation of hydrostatic equilibrium to the ICM, the temper- 
ature and density profiles are obtained once a single value is 
known, or selected. Thus, we fix the pressure in the outer- 
most bin to the value that allows the resulting deprojected 
temperature profile to match the observationally determined 
cluster temperature of 6.3 keV. This comes from the spectral 
analysis (a Raymond-Smith model plus a cooling fiow com- 
ponent) of the data collected from the Large Area Counter 
on the GINGA satellite (Allen et al. 1992). The gravitational 
potential functional form is defined by either of two difi'erent 
models for the dark matter distribution: (i) the true isother- 



mal sphere (Binney & Tremaine 1987); (ii) the Navarro- 
Frenk-White dark matter profile [pdm = Ps2;~^(l + 
Navarro et al. 1995] . A comparison between them for differ- 
ent input parameters is shown in Fig. ^. There, the [a, Tc] 
values for the case of the NFW potential are the velocity dis- 
persion of the dark matter under the isothermal condition 
1^, and the scale radius, as defined in Navarro et al. (1995). 

A gravitational contribution from the cD galaxy 
NGC 1275 has also been included using a de Vaucouleurs 
(1948) model with a velocity dispersion of 400 km s~^ (from 
the optical L — a'^ relation on the brightest cluster members 
in Malumuth & Kirshner 1985) and an effective radius of 21 
kpc (Schombert 1987). 

In Fig. ^, as an example of the information obtained by 
deprojecting a surface brightness profile, we show the results 
from the off-axis observation of the North of the Perseus 
cluster. 

The standard deviation confidence limits on the results 
summarized in Table 3 are estimated by perturbing the sur- 
face brightness profile ICQ times, according with the Poisso- 
nian error on the counts in each radial bin. For the cooling 
radius, rcooi, and the mass deposition rate, M, we quote the 
median values, with the respective 10th and 90th percentile 
limits, that come from the 100 perturbations. 

The gas mass, Mgas, and the total gravitating mass, 
Mtot , are obtained by integrating the gas and the dark mat- 
ter density, respectively, over the volume. The gas fraction, 
/gas, is defined as the ratio between the total gas mass and 
the gravitating mass within the radius 7?out, where the sur- 
face brightness is larger than its error by, at least, a factor 
of 3. 

The two dark matter profiles that we adopt, the true 
isothermal sphere and the Navarro- Frenk- White profile, give 
good agreement in describing the physical condition of the 
ICM under the hydrostatic and isothermal assumptions, 
apart from a systematic reduction (about 10 per cent in 
the averaged values) in the mass deposition rate when using 
a NFW dark matter profile rather than the BT profile. 

However, given the low resolution that we adopt, to 
avoid the technical limits on the observations (i.e., masked 
ribs, off-axis PSF) and to ensure a good signal-to-noise ratio 
in the outer part of the cluster, we do not investigate further 
the cooling flow region, apart from emphasizing that the 
results are roughly consistent with each other, and also with 
recent work that indicates a cooling radius of 0.2 Mpc and 
Mof about 500 Mq yr"^ (cf. Allen & Fabian 1997, Peres et 
al. 1998). 



2.2 The colour analysis 

Given the low spectral resolution of the PSPC [AE/E = 
0.43(_E;/0.93keV)-° '^; ROSATXJseis' Handbook 1994 []], we 
deflne three independent energy bands (or colours) where 



t = GAf20o/(2r2oo) = 50H^r^QQ, where M200 , r2oo , and G 
are the dark matter mass at the over-density of 200, the radius at 
which this over-density is reached, the Hubble constant and the 
Gravitational constant, respectively 

t also available at the URL 

heasarc . gsfc . nasa.gov /docs/ rosat /ruh /handbook /handbook, html 
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Table 3. Results of the deprojection analysis. The o-j^jj^g values come from the assumption Papcc = Amag.com where /3spec is the value 
obtained from spectroscopy and ftmag.corr is the best-fit result when a /3— model is adopted, as described in subsection 2.3.2. The errors 
on the gas fraction, /gas, consider the uncertainty on the temperature, kT, added in quadrature. Each column contains the two values 
obtained under the assumption that the cluster gravitational potential is described by either a true isothermal sphere (BT) or using the 
Navarro- Frcnk- White form (NFW). 

Region rc crdpr (o-imag) -Pout/10* Rout r-cool M /gas 

Mpc km s^i K cm-3 Mpc Mpc Mq yr"! 

BT, NFW BT, NFW BT, NFW BT, NFW BT, NFW 

North 0.25,0.35 860,860 1.7 1.93 0.20^2 ^2 > 0-20loj2 499tn. 474^11 0.31 ± 0.01, 0.31 ± 0.01 

West 0.45,0.60 830,860 1.6 2.29 0.21^0 03, 0.201°;^° 554^51,487^30 0.37 ± 0.01, 0.40 ± 0.01 



(557^26) 



South 0.60,1.80 970,1150 2.1 1.69 0.20+° j2, 0.20+o J° 489+45,452+2^ 0.21 ± 0.01, 0.22 ± 0.01 

East 0.50,0.70 960,1000 1.7 1.93 0.22tnn!!, 0.20lnn° 6291", 543^^3 0.24 ± 0.01, 0.26 ± 0.01 



{1029tll) 



Centre^off 0.40,0.70 890,960 1.2 2.29 0.21+°?,!^, 0.20tnno 519^?,, 477^1 0.30 ± 0.01, 0.29 ± 0.01 



(987^11) 



the analysis will be performed: (i) a soft band, which con- 
sider photons with energy enclosed between 0.52 and 0.90 
keV (formed by the sum of the R4 and R5 ROSAT bands; 
cf. table 1 in Snowden et al. 1994); (ii) a medium band, with 
an energy range of 0.91-1.31 keV (R6 T^OS'^T band), which 
contains most of the Fe-L emission produced by gas at about 
1 keV; (iii) a hard band (corresponding to R7) between 1.31 
and 2.01 keV. These 3 bands from on-axis PSPC response 
curve are shown in Fig. ^ 

We make use of the ratios among these colours to probe 
the physical properties of the Perseus cluster, once the ob- 
served ratio values are compared with the behaviour of the 
same ratio in a simulated cluster spectra. This kind of anal- 
ysis using the ROSATVSPC colours is also comprehensively 
discussed in Allen & Fabian (1997) in a study of the cooling 
flows properties of a sample of clusters. 

Predicted colour ratios were calculated in XSPEC (vers. 
10; Arnaud 1996) using as reference model the plasma code 
of Mewe, Kaastra and Liedahl {MEKAL model; Kaastra 
1992, Liedahl et al. 1995), with a metallicity, Z, of either 
40 or 70 per cent the Solar abundance Zq, and subjected 
to the photoelectric absorption provided from the model of 
Morrison & McCammon (1983). We present in Fig. |ll|, for 
a wide range of ICM temperature and Galactic absorption, 
the theoretical ratios obtained fixing one of the parameters. 

Firstly, these plots show that it appears difficult to de- 
tect any gradient in metallicity [note that the squares and 
the diamonds overlap]. They also show a weak dependence 
of the ratios upon temperatures higher than 3 keV. More 
meaningful are the steepening of the ratios at low temper- 
atures, particularly in the medium/hard [(91-131 PI)/(132- 
201 PI)] colors ratio, and the decrease in soft/hard [(52-90 
PI)/(132-201 PI)] when a higher column density is consid- 
ered (cf. the panels with triangles and diamonds). We also 



note that, if A^'h is assumed much lower than its true value, 
this effect hides the steep increase due to the presence of a 
lower gas temperature, particularly in the ratio where the 
softest band is involved. 

We performed a colour analysis along the four direc- 
tions (North, West, South, East) where the cluster has been 
observed. All the four offset observations allow us to extract 
the surface brightness profile, S{r), out to 80 arcmin from 
the X-ray centre. We use bins of [4,5,6,...] arcmin to increase 
the signal-to-noise ratio at larger radii (cf. Fig. [l2| . 

The background was selected in a 45°— sector opposite 
the peak of the X-ray emission and calculated as an average 
of the counts s~^ arcmin"'^ present between 40' and 45'. A 
Poissonian error is applied to the counts of photons present 
in the selected region. 

The errors, tcoi, on the colours ratio, ^coi, come from 
the propagation, in quadrature, of the l-cr uncertainty both 
in each radial bin and in the background: 

^ ^ ^ -S*! - 5'bkg,l £col ^ '^i + '^bkg.i 

S2 — 5'bkg,2 ' S^ol {Si — Sbkg,i)^ ' 

i — 1,2 

where S and e are the surface brightness value and the re- 
lated error, respectively, the index 'bkg' denotes the back- 
ground value. The sum is performed on the 2 energy bands 
considered. These errors increase towards the outer regions, 
becoming dominant at radius larger than 50 arcmin (around 
60 arcmin the signal to noise ratio, S/N, becomes ^ 3 for 
bin size of [4,5,6,...] arcmin; cf. Fig. [13 for the dependence 
of S/N on the bin size.). 

We describe below the results of the colour analysis per- 
formed through the ratios Scoi{r). Here we note that any ap- 
parent gradient has to be considered with respect to the er- 
ror bars. We observe a slight dependence on the background 
of the colour profile in the outer region, whilst the a devi- 
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Figure 9. Plots of the deprojection results of the northward observation of the Perseus cluster, when the potential for a true isothermal 
sphere is applied. The dotted lines in 'Temperature' and 'Cooling Time' plots represent the adopted Tgas of 6.3 keV and the Hubble 
time, respectively. In Table 3, we list the input and output values of the deprojection analysis. 



ations (showed in lower panels of Fig. |l3| - |l6|) are consistent 
whichever background area is selected. 

Our results are also compared to the conclusions from 
the (2.5-10 keV)/(l-2.5 keV) colour analysis on the Spartan 
1 data (Snyder et al. 1990). 

North The Spartan 1 results show in this direction a 
slight increase in the (2.5-10 keV)/(l-2.5 keV) ratio mov- 
ing away from the cluster centre, that they claim is related 
to the observed higher Galactic absorption to the NE of 
Perseus (cf. Fig. ^ ) . 



In our analysis (cf. Fig. |l^), we observe a 3(t fluctuation in 
the core region due to the cooling component. Apart from 
other 2(7 fluctuations around 20 arcmin, a plasma model 
with kT = 6.3 keV and iVn = 15.2 x 10^° atoms cm"^ is a 
good description of the data. 

West In this region. Spartan 1 does not show any colour 
gradient. 

A cool component (;$ 3 keV) is present at high significance 
in the core (panel b). The negative 3a deviation in the sec- 
ond bin can be explained only with an increase by a factor of 
about 2 in the absorption with respect to the value quoted 
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Figure 11. The diamonds connected by solid line represent the (T = 6.3 keV; AfH,20 = 14.9; Z/Zq = 0.4) case for the two different 
ratios considered, when either the temperature is varied between 2 and 7 keV (upper panels) or AfH,20 is in the range 10—30 (lower 
panels). In the panels at the top, the squares show the (varies; 14.9; 0.7) choice. In the panels at the bottom, the triangles represent 
the (3.0; varies; 0.4) case. These plots show that it appears difficult to detect any gradient in metallicity [note as the squares and the 
diamonds overlap]. They also show a weak dependence of the ratios upon temperatures higher than 3 keV, and a significant decrease 
when an higher column density is considered (cf. the panels with triangles and diamonds). 



of iVH,20 = 12.6 (Fig. 0). This 3a deviation is still evident 
between 15 and 30 arcmin, where the absorption remains 
high (cf. also Fig. 0) and GINGA (Arnaud et al. 1994) de- 
tects clumps at temperatures of about 9 keV. Both of these 
effects decrease the colour ratios, even if the excess in ab- 
sorption is dominant and requires roughly a AA'^h,20 ~ 3. 
Moreover, between 35 and 50 arcmin, there is a drastic in- 
crease in the ratios due to the presence of a cool component, 
corresponding to the region around the galaxy IC310. 

South Snyder et al. (1990) observe a very flat colour pro- 
file. 

Again, the colour analysis underlines the cooling region 
within 4 arcmin and an excess in absorption within 10 ar- 
cmin of about 4 X 10^" cm~'^ with respect to the quoted value 
of A'h,20 = 14.3. Furthermore, we obtain consistency with a 
kT = 6.3 keV above 10 arcmin, from panel (b) (Fig. |l5[). 



East At radii greater than 25 arcmin , Snyder et al. (1990; 
cf. their fig. lid) observe a softening of the cluster spectrum. 
With an estimated temperature of 2.5 keV, the core presents 
the larger deviation from the assumed isothermal model. 
Then, at least in the first 3 bins, we are overestimating A''h 
[from the temperature in panel (b), we require something 
more close to A'h,20 ~ 14]. We also note that the colour ratio 
behaviour is more consistent with a cooler ICM (Fig. |l6| 
shows less dispersion for the solid line that represents the 
gas at 3 keV) between 20 and 50 arcmin. Here, a slight 
positive gradient in the colour profiles is present between the 
third and sixth bin (i.e. 9-39 arcmin), underlying a softening 
of the emission where the excess in surface brightness is 
detected (cf. also Fig. |^ and 

To summarize, we observe variations in X-ray colour 
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Figure 13. The colour ratios for the North region. The Galactic absorption value used is quoted in Table 1. (Upper panels) The colour 
ratios 5coi('') are compared with the predicted value at 3 keV (solid line) and at 6.3 keV (dashed line). (Lower panels) The deviation (in 
(t) from the predicted values (3 keV: solid line; 6.3 keV: dashed line) are plotted. In these panels, the dotted line shows the zero level 
and the dash-dot-dot line the ±3(t level. 
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Figure 14. As Fig. ha for the West region. 
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Figure 15. As Fig. ^ for the South region. 
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Figure 16. As Fig. hsifor the East region. 
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Figure 10. These are the Effective Areas of the 3 energy bands 
considered. They are extracted from the Redistribution Matrix 
File pspcb_gain2_256.rmf . 
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Figure 12. Signal to noise ratio for Soft/Hard band ratio in the 
East with bins size of 30 arcsec, 90 arcsec and 180 arcsec (upwards 
solid lines). The dashed line is for bins size of [4,5,6,...] arcmin 
and its maximum is about 90 in the first 2 bins. The dotted line 
represents the threshold of SjN = 3. 



across the four regfons. There is a clear temperature de- 
crease in the core associated with the cooling flow. At larger 
radii from the centre the changes are best explained by small 

variations in A^h. There is no evidence for an overall system- 
atic temperature gradient up to 50 arcmin (~ 1.5 Mpc). 



2.3 The /3 problem 

Cavaliere & Pusco-Femiano (1976) pointed out that the scale 
height for galaxies and gas in a cluster may be calculated 



Table 4. Best /3-model fit values from previous work. The er- 
rors are at the 90 per cent confidence level. The labels in col- 
umn re/ indicate: Cr97 (Cruddace et al. 1997), Sc92 (Schwarz 
et al. 1992), Sn90 (Snyder et al. 1990), JF84 (Jones & Forman 

1984), Br81 (Branduardi-Raymont et al. 1981); Cr97 and Sn90 
apply two different method on the same data from the Spartan 
I: Cr97 bin the spatial and the spectral counts distribution and 
fit a three-dimensional model of the cluster structure, obtained 
by the Raymond-Smith code of the plasma emissivity; Sn90 use a 
/3-model on the surface brightness profile. They quote a value 
of 287.2 (with 273 degree-of-freedom) and 20.1 (19), respectively. 
The latter method is also used by Sc92 on the ROSAT All Sky 
Survey data, and by JF84 and Br81 on the Einstein IPC data (cf. 
also fig. 8 in Sn90). In Sc92, both errors and values are not 
quoted. In JF84, the value is not indicated. Also in Br81, no 
X^ value is indicated and the core radius is fixed to the one deter- 
mined from the optical galaxy distribution [using a King model, 
Kent & Sargent (1983) find a core radius of 11 arcmin]. 



Region 




rc 


ref 


r < 50' 

10' < r < 78' 


'•-'-0.09 

0.52 


13.3'_+,^-«' 
1.9' 


Cr97 

Sc92 


8' < r < 50' 
3.2' < r- < 32' 


r. -9-H0.08 
^■'^-0.05 

0.57 ±0.03 


8.9' ± 1.7' 


Sn90 
JF84 


10' < r- < 17' 


0.59 ±0.02 


8' (fix) 


Br81 



from optical velocity dispersion Uopt and X-ray spectral tem- 
perature kT as 



/^spec — 



(2) 



and, for an isothermal spherically-symmetric X-ray emis- 
sion, from the surface brightness profile in the functional 
form (when r 3> rc) 

0.5 3/3iniag 



S{r) = So 



1 + 



(3) 



Edge & Stewart (1991) calculated from the EXOSAT 
data on the Perseus cluster a /Sspec of 1.78^0 31 . Using 
the more recent reference values (with corresponding l-cr 
errors) on tTopt of 1026lgf km s~^ (Fadda ct al. 1996) 
and kT = 6.33to\ll keV (Allen et al. 1992), we obtain 
/Sspec = 1.091q ;j4. The significant decrease comes mainly 
from a much more careful determination of ffopt , once sub- 
structures, gaps in galaxy velocity distribution and more 
robust statistical tools are considered. 

However, this estimate is much larger, for example, than 
Amag of about 0.52, obtained from the azimuthally aver- 
aged profile of the cluster surface omission (Schwarz ct al. 
1992; but also cf. Table 4). The disagreement between flspoc 
and /3imag values is the reason for the so-called ,5 problem 
in the Perseus cluster (cf. also Allen et al. 1992). In partic- 
ular, the non-symmetric emission from the cluster (as first 
pointed out by Branduardi-Raymont et al. 1981) traces a 
recent merger, that, together with the excess of gas to the 
East and of galaxies to the West, seems to emphasize the 
dependence of the /9— problem on the sector analysed. In 
fact, Schwarz et al. (1992) observed a value of Amag ~ 0.5 
to West and North, and of 1.27l|] js East (in their paper, 
there is no mention of the (3 value to South). 

We re-examine this issue, using our offset profiles out 
to about 80 arcmin from the X-ray peak. They are binned 
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Figure 19. The 68, 90 and 99 per cent confidence contours in 
the rc — /3 plane for the /3-model applied on the 0.5-2 keV surface 
brightness profile selected between 10 and 80 arcmin (cf. the last 
point on the thickest solid line in Fig. ha). 



in width of [ 4, 3, 2, 1, (all the remaining at) 1] arcmin, due 
to the PSF behaviour in the outer part of PSPC. 

First, an attempt is made to fit the profiles over the 
whole range of [ 0, 80] arcmin (cf. Table 5). Here, we note 
that, applying a model, the goodness-of-fit is generally 
very poor, in a sense. This is due to the fact that the sta- 
tistical Poissonian errors related to the total counts present 
in each radial bin are particularly small (cf. Fig. They 
do not consider any fluctuations present in each bin due to 
azimuthal variations within each sector. The contribution 
from these azimuthal fiuctuations will be discussed in the 
next section. 

Furthermore, as just noted in previous work (cf. Ta- 
ble 4), there is enhanced emission above the central part in 
all models. Thus, we have also analysed the surface bright- 
ness profiles avoiding the central 10 arcmin affected by the 
presence of the large cooling fiow (as shown in the depro- 
jection analysis, it has a typical rcooi '~ 0.2 Mpc, equivalent 
to about 6.5 arcmin; cf. also the colour ratios). As indicated 
in Table 5, the reduced x'^ value decreases significantly but 
still remains too large to be statistically acceptable. 



2.3.1 Estimated error in the surface brightness 
mean value 

Previously, we have measured the error on the surface 
brightness as the Poissonian error related to the total counts 
in each radial bin. This is a statistical error that does not 
consider any fiuctuation on the mean value (i.e., the surface 
brightness value) due to azimuthal variations intrinsic to the 
cluster. 

The contribution expected from these fiuctuations is 
particularly high in a non-relaxed cluster such as Perseus. 



Figure 20. This is a 45'-radius map of the residuals significant 
at the 95 per cent confidence limit in terms of the probability 
distribution as drawn by equation (5). It is obtained after the 
mosaic of the four exposure-corrected image has been subtracted 
from the cluster emission simulated by a /3-model with param- 
eters Tc = 17.0' (68 per cent uncertainties range: 13.5, 21.3), 
/3 = 0.81 (0.75, 0.90), and 5o = 2.31 x 10"^ counts s-^ arcmin'^. 
With a reduced of 0.05, these values provide the best-fit to the 
Centre-ofi' surface brightness profile, between 10 and 80 arcmin 
(cf. Table 5). 



We have estimated the contribution from the dispersion of 
the mean values (assumed normally distributed) of the sur- 
face brightness calculated in A^sect azimuthal sectors. This 
is done in each off-axis image (where the maximum angu- 
lar size is about 90° at the maximum width). Choosing 
Ascct = 4 (i.e. angular sector of about 15°), we obtain a 
relative error of about 15 per cent over the radial range 
[ 0, 80] arcmin, i.e. about 10 times larger than the simple 
Poissonian error. The surface brightness profiles with the az- 
imuthal fiuctuations added in quadrature to the Poissonian 
errors are shown in Fig. |l^. 

After fitting these surface brightness profiles, we obtain 
values of reduced x'^ which are low due to the large, but real, 
errors now present in the data (cf. Tab. 4). 

Of course, the value of is reduced by increasing Asoct, 
as fiuctuations further dominate. On the other hand, it ap- 
proaches a more statistically meaningful value of about 1 
when both the angular region considered has width ^ 15° 
and only Poissonian error is considered. Then, the contri- 
bution from fiuctuations is still larger by a factor of 3 and 
more. 

In the following analysis, we quote the best-fit results 
obtained over each 90°-offset region (cf. Tab. 4), that can be 
considered as representative of the variations along different 
directions. 

In Fig. |l8[ we plot the measured /3imag and core radius 
for each sector, at different outer radii and for each of the 
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Table 5. Best fit values of the surface brightness profiles and the corresponding reduced x'^i x2i adopting difi^erent models. The values 
quoted are (scale radius -in arcminutes-, slope). Particularly, we present (re, /9) for the /3-model (rcut; ''p, 7; fc 13) for the 'power law 
+ /3' model (pow + /3); ('"cut ; ''pii 7i ; ''P2 1 72) for a two broken power laws (2 pow); (rs,r)) for the NFW gas density profile (NFW gas); 
(rv, ce) for the de Vaucouleurs law (de Vauc). 







Region 


model 


values 










only Poissonian error 




North: 


10' 


< r < 80' 




( 8.8, 0.65) 


1.90 




0' 


< r < 80' 


/3 


( 3.6, 0.58) 


6.67 


West: 


10' 


< r < 80' 




(12.5, 0.60) 


12.87 




0' 


< r < 80' 


/3 


( 2.0, 0.51) 


51.15 


South: 


10' 


< r < 80' 




(25.3, 1.05) 


2.60 




0' 


< r < 80' 


/3 


( 2.2, 0.53) 


71.12 


East: 


10' 


< r < 80' 




(24.8, 1.23) 


8.55 




0' 


< r < 80' 


/3 


( 4.4, 0.60) 


70.89 


Centre-off: 


10' 


< r < 80' 




(17.5, 0.85) 


15.42 




0' 


< r < 80' 


/3 


( 3.4, 0.57) 


122.38 






Poissonian 


error and azimuthal fluctuations 




North: 


10' 


< r < 80' 




(9.5,0.67) 


0.12 




0' 


< r < 80' 


13 


(4.3,0.62) 


0.17 








pow + (3 


(10.0; 3.5,1.21; 9.5,0.67) 


0.14 








2 pow^ 


n ^ D- 9 1 9S- 4 1 

^^J-O.U, ^.U,-L.^O, ^.0,-l-.0( J 


0.18 








NFW gas 


(15.6,8.87) 


0.22 








de Vauc 


(0.017,0.29) 


0.13 


West: 


10' 


< r < 80' 


13 


(12.9, 0.61) 


0.37 




0' 


< r < 80' 


[3 


(2.7, 0.52) 


1.00 








pow + (3 


(27.0; 3.3,1.38; 67.6,1.67) 


0.11 








2 pow 


(26.0; 2.0,1.38; 4.3,1.77) 


0.34 








NFW gas 


(25.3,7.84) 


1.08 








de Vauc 


(0.001,0.23) 


0.54 


South: 


10' 


< r < 80' 


13 


(25.8,1.10) 


0.19 




0' 


< r < 80' 


[3 


(21.7,0.99) 


0.72 








pow + (3 


(12.0; 3.4,1.39; 25.6,1.07) 


0.23 








2 pow 


(26.0; 2.0,1.33; 6.1,2.39) 


0.32 








NFW gas 


(60.7,12.79) 


0.66 








de Vauc 


( 3.33,0.67) 


0.54 


East: 


1 r\f 
W 


< y < 80' 


13 


(14.5, 0.88) 


U.4U 




0' 


< r < 80' 


P 


(10.1, 0.79) 


0.64 








pow + (3 


(18.0; 3.5,1.22; 6.4,0.79) 


0.33 








2 pow 


(19.0; 2.0,1.22; 5.7,2.30) 


0.33 








NFW gas 


(28.5,10.82) 


0.52 








de Vauc 


(0.42,0.46) 


0.50 


Centre-off: 


10' 


< r < 80' 


13 


(17.0,0.81) 


0.05 




0' 


< r < 80' 


13 


(4.7,0.63) 


0.54 








pow + /3 


(10.0; 3.4,1.24; 17.0,0.81) 


0.09 








2 pow 


(22.0; 2.0,1.31; 5.2,2.12) 


0.14 








NFW gas 


(20.3,9.37) 


0.58 








de Vauc 


(0.045,0.33) 


0.13 



three energy bands. Thus, we check also the dependence of 
/3imag upon the colours, that we expect if any merging action 
is produced by substructure with a significant difference in 
temperature from the ICM. 

Using the /3— model for the analysis of the different sec- 
tors, we find a value around 0.65 in the North, and an in- 
crease from the exterior to the centre, from 0.45 to 0.65, in 
the West. 

The eastern and southern regions provide values of 
about 1 over all radial ranges considered. These values, 
as shown by Edge & Stewart (1991), indicate an on-going 
merger that seems to involve the eastern region (cf. again 
the conclusion in Schwarz et al. 1992) as well as the south- 
ern one. In particular, the soft profile steepens significantly 



around 40 arcmin to the East. The colour analysis shows 
a gradient that agrees with a softening of the emission (cf. 
Fig.|l|). 

Due to the correlation in the fi — Tc plane (cf. Fig. |l^ , 
the core radii emphasize a peculiar gas distribution in these 2 
sectors: they are about twice the ones estimated in the North 
and West. The presence of this correlation can also partially 
explain the disagreement between the parameters fitted in 
the previous work, where the lower core radii correspond to 
lower P (cf. Table 4). 

Fig. ^ also shows that, although there is a trend for 
the gas distribution to be fiatter in the medium band, the 
surface brightness selected in the energy range 0.9-1.3 keV 
shows a steeper profile in the East sector, underlining the 
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Figure 17. The [ 4, 3, 2, 1, (all the remaining) 1] arcmin plots with the errors obtained considering also the azimuthal fluctuations 
within each radial bin. These plots are fitted with the /3-model and the other models discussed in section 2.4. The results of the fitting 
procedure are presented in Table 5. 



ICM interaction with a clump of gas at temperature greater 
than 1 keV. 

The 0.5-2.0 keV band indicates a /3imag value consis- 
tent with the average of the 3 other values calculated in the 
different bands. 

The errors on the two interesting parameters, /Jimag and 
Tc, are shown in Fig. |l9|as contours plotted at 68, 90 and 99 
per cent of confidence level (i.e. x^-value regions of Xmin + 
2.30, Xmin+4.61 and Xmin+9-21, respectively; e.g. Press et al. 
1992, page 687). They are estimated with the fit performed 
on the surface brightness profile selected between 10 and 
80 arcmin in 0.5-2.0 keV band. No agreement is evident, at 
the 99 per cent confidence level, among the different sectors 
considered. Thus, the complexity of the area in terms of 
the gas distribution is once more emphasized. However, a 
good representation of the region is provided through the 
azimuthally averaged Centre-ojf profile. 

With these new constraints on the /Jimag values in the 
different sectors of A426, we can address the issue of the 
/3-problem. 



2.3.2 Resolution of the P problem 

Following Bahcall & Lubin (1994), and the recent results on 
N-body simulations of clusters of galaxies (cf. Thomas et al. 
1998), we first consider a correction to the /3imag value. The 
/3-model, as described by equation (3), has been introduced 
assuming that the gas and the galaxies are in hydrostatic 
equilibrium with the cluster potential, and that the galaxy 
distribution is well described by the King approximation 
(1962) to the isothermal sphere. This approximation is pro- 
portional to at the outer radii, and does not match the 
dark matter dependence of r"'^ '*, which is determined from 
N-body simulations (cf. Fig. ^). Then, the hypothesis that 
the optical light traces the mass, i.e. that the galaxy number 
density profile describes the gravitational profile of the clus- 
ter over the virialized region (as the recent CNOC results 
also seem to confirm; cf. Carlberg et al. 1997), indicates a 
way to correct the /3imag value to give the correct dependence 
on the galaxy distribution of the radius: 

Pgas(r) = Pgai(r)'''"''s oc r"^'''»'"5 cc r-^-'^'f'^s.-"" , (4) 



i.e., /3imag,corr — 1.25A n 

© 0000 RAS, MNRAS 000, 000-000 



16 



North 



0.8 - 



0.6 - 




North 




West West 




30 40 so 60 70 80 90 30 40 SO 60 70 80 90 





Figure 18. These plots show the results on the P value and the core radius, Tc, when the /3— model is fitted on the profiles binned in [ 
4, 3, 2, 1, (all the remaining) 1] arcmin and selected between 10 arcmin and different outer radius in the 3 energy bands: soft {triangles 
joined by a dotted line); medium (dashed line and squares); hard (solid line and diamonds). The thickest solid line represents the result 
when the /3— model is applied to the 0.5-2.0 keV exposure-corrected image. At the last point on this line are related the error-bars drawn. 
These plots show a correlation between /3 and rc (see also Jones & Forman 1984). 
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Thus, we have compared /3spec ~ 1.091q ^4 with the cor- 
rected values of /Jimag; 0.84, 0.76, 1.37, 1.10, 1.01 (with 68 per 
cent errors of about 10 per cent) for the North, West, South, 
East and Centre-off, respectively. There is now a reasonable 
agreement between /3imag,corr and Pspcc, which eliminates the 
P problem in the Perseus cluster. 

Furthermore, assuming the equation /3spec ~ Amag.corr 
is valid, we can calculate a velocity dispersion, o"imag, us- 
ing /3imag values (from the 0.5-2.0 keV band selected in the 
radial range 10-80 arcmin) and an isothermal temperature 
of 6.3 keV. These values are quoted in square brackets in 
Table 3, and show a good agreement with the deprojection 
estimates, (Jdpr. AH the velocity dispersions calculated so, 
are also marginally consistent with the optically determined 
velocity dispersion, (Topt, of 1026l^f km s-\ 

Using the best-fit results from the /3— model, we inves- 
tigate in Fig. the residuals (significant at the 95 per cent 
confidence level) for a cluster simulated by using a /3 - model, 
S{x,y), compared to the raw data, I(x,y). Both S{x,y) and 
I{x,y) are in counts pixel"^. Then, the significance is as- 
sessed in terms of a distribution: 



2 

X = 



N^XNy 

3 = 1 



S{xj,yj) 



(5) 



where Nx and Ny are the numbers of pixels in the rows and 
columns in the selected sectors of 22? 5 in annular circular 
bins of [ 4, 3, 2, 1, (all the remaining) 1] arcmin. 

From this map, apart from the central cooling region 
which contributes to the total luminosity Lboi = 2.8 x 10*^ 
erg s~^ about 50 per cent (cf. the deprojection results in 
in Fig. the eastern excess extends to 40 arcmin from the 
centre with an estimated luminosity in the residuals of about 
8 X 10*'^ erg s~^. From the L — T relationship in the groups 
and clusters sample of Mulchaey and Zabludoff (1998), this 
luminosity corresponds to about 2 keV, as we also conclude 
from the colour analysis. Other extended emission appears 
related to the galaxy IC310, at about 40 arcmin westwards 
(IC310 is the masked point source; see also Cruddace et al. 
1997). This clump has an estimated luminosity of ~ 1 x lO**^ 
erg s~^ with a corresponding temperature of 1 keV. This is 
further support for the drastic increase in the colour ratio 
values for the same region is attributed to an excess of cool 
gas, as seen in Fig. ^ Both of these clumps appear as groups 
still merging with the cluster. 

Other poorly defined discrepancies with an azimuthally 
averaged /3-model profile appear in the South. They can be 
part of the eastern excess and are responsible of the /3 value 
of about unity measured in that sector. 

We stress that these residuals depends on the assump- 
tion of the model. As shown in Section 2, the residuals are 
not significant if the elliptical isophotes are allowed to have 
free centroids. Furthermore, modelling this elliptical isopho- 
tal intensity profile with a /3-model in the range [10-80] ar- 
cmin, we obtain P = 0.63 ± 0.01 and Vc = 9.3 ± 0.2 arcmin, 
for a reduced of 33.6. These values agree well with those 
observed for the sector to the North, that appears as the 
most relaxed region in the cluster. 



2.4 On the fitting of the surface brightness profile 

The inadequacy of the /3-model in reproducing the gas profile 
over the whole radial range, in particular in the central part, 
has driven us to look for other possible functional forms. 
We consider the following models: 

(i) a power law plus a /3 model, to describe the inner 
cooling region and the outer part, respectively, with a radius 
of cut-off, rcut, as a free parameter: 



1-27 



if r < Vcu 

So ifr->rcu 



(6) 



where the exponent for the 'power law' model comes from 
the consideration that 



Pgas = PO 



Sb ~ So 



1-27 



(7) 



(ii) two density power laws, boken at radius rcut; 

(iii) a generalized form of the de Vaucouleur's law (1948): 



Sh = So exp 



(8) 



We also introduce a gas density profile obtained from 
the NFW dark matter profile when the hydrostatic equation 
under the isothermal condition is applied: 



Pgas = Po(l + X) 



(9) 



where x — r/rs and r) = A-KGpsr'i^mp/{kT) (cf. Makino, 
Sasaki & Suto, 1998, and Ettori & Fabian, 1998). The sur- 
face brightness profile is then obtained by numerical inte- 
gration of the gas density profile in equation (9). 

We use the reduced provided by each model to as- 
sess their suitability. The best-fit parameters are listed in 
Table 5. From this table, the 'power law -I- (5 model' ap- 
pears as the more appropriate description of the data over 
all the range of radii. In terms of the _F-test, when we com- 
pare the 13 model with the 'power law -I- /3' fit, we obtain 
that it is still a good description of the northern data, but it 
becomes worst at 99 per cent of confidence level in the cases 
of the West, South, East and Centre-off. Once again (cf. the 
colour analysis and Fig. 18), the North appears as the more 
relaxed region in the Perseus cluster. 

We also note that the 'power law -I- /? model' highlights 
the merger of the clump of gas around IC310 in the outer 
region to West. In fact, the P value increases there to about 
1.4, compared to a value of 0.6 from the single /? model. 
Moreover, the cut-off radius is 26 arcmin to West and 18 
arcmin to East (around 12 in the other sectors), underlining 
the elliptical configuration discussed in Section 2. 

However, the /3 model is a poor description when it is 
also compared to the models (ii) and (iii), but it is better 
than the gas profile obtained from the NFW law. In fact, 
this latter profile shows the greatest disagreement with the 
data. We check elsewhere on the consistency of the gas dis- 
tribution obtained from NFW dark matter potential with a 
larger set of X-ray cluster profiles. 

The slope in pgas oc r~" is in the range [1.22, 1.38] in the 
core region and become steeper to [1.77, 2.39] in the outer 
part of the Perseus cluster. There, the slopes from the two 
broken power laws are consistent, within an average error of 
15 per cent, with the asymptotic behaviour of the /3 model 
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Figure 22. The dark matter density profiles as obtained from 
the best-fit values of the NFW gas profile are the dashed line 
{Centre- ojf) and the dotted lines (the other four sectors). The 
crosses represent the result from the (3 model, the triangles are 
from the deprojection analysis, and the squares what we obtain 
applying two broken power laws (we plot just the external values 
for a more clarity) . All of them are obtained in the case of Centre- 
off. 



(apart from the southern region where the discrepancy is 28 
per cent). 



3 THE CONSTRAINTS ON THE MASS 
PROFILE 

The form of the gravitational potential remains the most 
uncertain element of the deprojection analysis, since the ob- 
servational constraints are not completely sufficient. How- 
ever, some conclusions on the shape of the dark matter pro- 
file can be achieved under the assumptions of hydrostatic 
equilibrium and isothermality of the gas distribution in the 
cluster. Then, the total gravitating mass, responsible for the 
external potential that maintains the gas mass within a ra- 
dius R in equilibrium, is written as 



47r /9dm r dr = 



dr 



(10) 



This equation, together with the perfect gas law, is used 
in Fig. ^ to describe how the total gravitating mass depends 
on the radius, applying both the deprojection technique and 
the /3-model. The results indicate consistency between the 
total mass, estimated by these methods, and previous esti- 
mates (when the errors are properly considered). In the same 
figure, the mass distribution, as obtained through the de- 
projection analysis, emphasizes the over-density in the east 
sector where the merging group acts. 

In Fig. we plot the dark matter density profile as ob- 
tained from (i) the best-fit (rj, rs) values of NFW gas profile. 



once the definition of rj is used with a plasma temperature of 
6.3 keV; (ii) the best-fit slopes 7 from the two broken power 
laws model, given that Mtot oc 7 kT r and pdm oc 7 kT r~^\ 
(ui) the /3-model, in which Mtot (x (5 kT /{r'^ -\- rl) and 
Pdm oc P kT {r'^ +'ir1) / {r'^ -\-r'^)'^ ; (iv) the result of the depro- 
jection analysis using the potential from the true isothermal 
sphere. Here we note that both the /3 and the power law 
models indicate a dark matter density profile that is flatter 
than N-body simulations profile, with a dependency upon 
radius of slope 2, instead of about 2.4. However, on this ra- 
dial range, this affects slightly the determination of the total 
mass as shown in Fig. ^ 

From fitting the core with a power law, it appears 
more convenient to use a steeper profile than the NFW one, 
whereas the King profile provides a plateau (cf. also Fig. 
Furthermore, the best-fit values define 7-200, the radius at 
which the average over-density in the cluster is 200, equal 
to about 2.7 Mpc. This means that we are mapping, in our 
analysis, the region where the over-density is about 300 or 
more. 

The results in Fig. ^ show good agreement within the 
region 0.4 - 0.8 r2oo, i.e. above the core region that de- 
pends upon the functional form of the potential adopted 
and within the radial limits of our observations. 

These constraints on the total mass are applied to de- 
scribe the distribution of the gas fraction, /gas, shown in 
Fig. ^ in dependence both of the central azimuthally aver- 
aged observations and of each sector considered. We obtain 
a gas contribution to the total mass of 30 per cent at 2.3 
Mpc. This value is consistent with the SpaHan 1 estimate 
at 1.5 Mpc of about 24 per cent (from Fig. 3 in Cruddace 
et al. 1997). 

Here we note that if any decrease in the average tem- 
perature of the cluster plasma is detected in the outskirts, 
as appears common in clusters from the recent analyses of 
ASCA data by Markevitch et al. (1997), then the gas frac- 
tion obtained assuming an isothermal profile should be in- 
creased significantly due to the corresponding reduction of 
gravitating mass. However, from the colour analysis, we ob- 
serve that the gas temperature in Perseus does not decrease 
1^ with respect to the reference value of 6.3 keV up to 50 
arcmin (about 10 rc), at least in the North and South. In 
the other two sectors, we detect cool emission from extended 
and identified components, probably due to a merger that 
could affect the hydrostatic assumption. 

The latter two excesses in emission are also seen in the 
plots of the gas fraction (Fig. where the distribution of 
each sector is compared to the azimuthally averaged value. 
The eastern excess, between 0.5 and 1 Mpc, corresponds 
to the merging group, while the excess in the West, where 
the gas fraction reaches a value of about 40 per cent at 2.3 
Mpc, shows a discrepancy with the azimuthally averaged 
value over all the radial range outside 1 Mpc. This is also 
shown in Fig. where the contribution from the extended 
emission around IC310 becomes significant. 



§ using the same analysis, it is difficult to put any upper limit; 
cf. Fig. 
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Figure 21. These plots present the mass distribution how has been obtained applying the deprojection technique, under the hydrostatic 
and isothermal assumptions, (left) The map is a combination of the four off-axis observations results on the total gravitating mass, (right) 

The mass profiles obtained by applying both the deprojection technique and the best-fit results from the /3-model on the azimuthally 
averaged surface brightness profiles are here compared with the previous estimates: about 5 X at R=1.3 Mpc {cross; Eyles et 

al. 1991); 6.7 x at 1.5 Mpc (square: from Cruddace et al. 1997 after we adopt the redshift of 0.0178). The last point at 2.3 Mpc 

reports a value of 1.19^Q'j^2 ^ 10^^ Mq, for the range of the 10th and 90th percentile on 100 Monte-Carlo simulated estimations (see 
Sect. 4). 



4 CONCLUSION 

Our main conclusions are as follows: 

(i) The distribution of the ICM in the Perseus cluster is 
complex, being more homogeneous in the northern region, 
with clear excess in the eastern part. We can not confirm 
any temperature higher than about 6 keV in the North, or 
North- West areas, as claimed from the ASCA results. 

The analysis of the PSPC colour ratios indicates an 
isothermal profile around 6 keV within 50 arcmin, a sig- 
nificant decrease in the cooling fiow region, and fiuctuations 
both in the eastern (between 20 and 50 arcmin from the 
X-ray center) and in the western areas (around IC310, 35- 
50 arcmin from the center), where wc detect cool extended 
emission with luminosity of few times 10**^ erg s~^. This re- 
sult is also supported from a /3-model analysis and shows 
the presence of groups still merging with the main body 
of the cluster. This merger event produces a disturbance of 
the gas within the gravitational potential, that can be mod- 
elled by elliptical isophotes with free centroids. This model 
then eliminates significant residuals from the regions of the 
clumps, emphasizing that these features must rather be en- 
tities within the cluster, and not foreground / background 
objects. 

Furthermore, we show that the observed discrepancy in 
surface brightness can not be justified simply in the basis 
of the gradient in absorption present over the region, even 
through some observed gradient is present as detected from 
our colour analysis. 

(ii) We use the /3-model to highlight the disturbed intra- 



cluster gas distribution in the East, South-East. The best- 
fit parameter /3imag is shown to be dependent on the outer 
radius in the extraction of the surface brightness profile, 
although it converges outside 50 arcmin (1.5 Mpc). The dis- 
agreement between /3spcc and ftmag disappears when /3spoc 
is properly evaluated with a recent estimate of the optical 
velocity dispersion that takes into account substructures in 
the velocity space, and /Jimag is carefully estimated in the 
outer part of the cluster with a proper correction for the 
slope of the dark matter density profile. 

Thus, we can claim that the /3-problem in the Perseus 
cluster is solved. 

(iii) The surface brightness profile in the Perseus cluster 
does not seem to be well described by a single /3-model over 

all the range [ 0, 80] arcmin. This is common for cooling 
flow clusters that also show merging in the outskirts, as in 
Perseus. 

After testing several different models, we find that the 
best statistical agreement between the data and the gas dis- 
tribution is obtained with a model with a power law fitted 
to the inner cooling region, and a /3 model for the outer 
component. 

We also apply of the gas profile obtained using the NFW 
dark matter potential in the hydrostatic equation. How- 
ever, this does not provide an acceptable description of the 
Perseus data, due to the no-relaxed state of the cluster. 

(iv) With the assumption of hydrostatic equilibrium, we 
place constraints on the non-luminous mass density profile 
from the best-fit results of the NFW gas profile, /3-model, 
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Figure 23. These plots show the gas fraction, /gas, as function of the radius R for the four sectors deprojected by assuming a Navarro- 
Prenk- White profile dark matter density profile (triangles joined by a dotted line). We also show, as comparison, the value for the gas 
fraction measured in the azimuthally averaged deprojected image, both using a BT profile (solid line) and a NFW profile under isothermal 
assumption (dashed line). (Cf. also Table 3 for the values of /gas quoted at Rout and the relative error). At 1.5 Mpc, Cruddace et al. 
(1997) calculate a value of about 0.24. 



two broken power laws model, and the conditions provided 
from the deprojoction analysis. 

All of these are consistent in defining r2oo at 2.7 Mpc. 

The total mass in cluster within 2.3 Mpc (0.85 r2oo) is 
1.19t^ la X 10^^ M0, where the errors are the 10th and 90th 
percentile calculated estimating 100 times the mass through 
the /3— model, once the temperature is randomly selected 
around its best value and the errors arc considered, the core 
radius is also randomly selected and f3 is fixed from the best 
fit on the surface brightness profile. The different models 
agree (within the respective error) on this estimate, even if 
a slight difference 20 per cent) is present in the slope 
of their dark matter profiles in the range of 0.4-0.8 r2oo, 
where our observations of the gas profile are available and 
not affected from the core region. 

(v) Using constraints on the gas and non-luminous mass 
distribution, we measure the gas fraction to be equal to 
about 30 per cent at 2.3 Mpc. This is similar to that ob- 
served in Coma and other clusters (White et al. 1993, White 
& Fabian 1995), and disagrees with the primordial nucle- 
osynthesis estimates of the baryon fraction in the Universe 



(e.g. Olive 1996), unless the ratio between the total density 
and the critical value, fio, matter, is less than 0.2. 
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